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ABSTRACT 

The nebula powered by the Vela pulsar is one of the best examples of an 
evolved pulsar wind nebula, allowing to access the particle injection history and 
the interaction with the supernova ejecta. We report on the INTEGRAL dis- 
covery of extended emission above 18 keV from the Vela nebula. The northern 
side has no known counterparts and it appears larger and more significant than 
the southern one, which is in turn partially coincident with the cocoon, the soft 
X-ray and TeV filament towards the centre of the remnant. We also present 
the spectrum of the Vela nebula in the 18-400 keV energy range as measured 
by IBIS/ISGRI and SPI onboard the INTEGRAL satellite. The apparent dis- 
crepancy between IBIS/ISGRI, SPI, and previous measurements is understood 
in terms of point spread function, supporting the hypothesis of a nebula more 
diffuse than previously thought. A break at ~25 keV is found in the spectrum 
within 6' from the pulsar after including the Suzaku XIS data. Interpreted as 
a cooling break, this points out that the inner nebula is composed by electrons 
injected in the last ~2000 years. Broad-band modeling also implies a magnetic 
field higher than 10 fiG in this region. Finally, we discuss the nature of the north- 
ern emission, which might be due to fresh particles injected after the passage of 
the reverse shock. 
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Subject headings: pulsars: general — pulsars: individual (PSR B0833-45) - 
ISM: supernova remnants — ISM: individual objects (Vela PWN) - - X-rays: 
general 



Introduction 



Pulsar Wind Nebulae (PWNe) are the non-thermal bubbles inflated by the winds 
of rotation-powered pulsars. Recent observations have allowed to arrange the variety of 
their morpholo gies in an evolutionary sequence resulting from the interaction with their 
surroundings (IGaensler fc Sland 120061 ) . A very complex phase occurs after the host rem- 
nant evolves into the Sedov- Taylor phase (~10 kyr after the pulsar birth), when a reverse 
shock propagates inwards in to the supernova ejecta and eventuall y collides with the PWN 
( Ivan der Swaluw et al.l 120011 ) . Located at a distance of 290 pc (IDodson et al.l 120031 ) . the 
PWN powered by the Vela pulsar (PSR B 0833-45, with spin- down luminosity E = 6.9 x 10 36 



erg s 1 , and characteristic age ~11.4 kyr; IDodson et al.ll2007l ) is the best example of a PWN 



in such a phase, which is thought to be typical of composite supernova remnants. 

The prime evidence of the pulsar wind is t he structured PWN reso lved by Chandra 
within ~1' fr om the pulsar (THerfand et al.l boOll: IPavlov et allboOli bpOflh . The X-ray flux 
measured by iMangano et al.l (120051 ) at varying integration radii shows that the Chandra 
PWN is surrounded by fainter emission up to at least 15'. At a larger scale, ROSAT found 
a X-ray filament (the so-calle d cocoon) extending ~45' south/south- west from the pulsar 
( iMarkwardt fc Ogelmanlll995l ). The cocoon also stands out amon g the radio filaments com- 
posing the large (~1.5° radius) Vela X nebula ( iRishbethl 119581 ). All these structures are 
embedded in the Vela supernova remnant, which has a radius of ~4°. 

The cocoon is ascribed to an asymmetric reverse s hock, which stripped the bulk of the 
particles from around the pulsar leaving a relic PWN (IBlondin et al.ll200ll ). This explana- 
tion has been confirme d by the detection of extended TeV emission matching the cocoon 
( lAharonian et al.ll2006l ). with the brightness peak offset f rom the pulsar , and of thermal 
X-ray emission suggesting mixing with the shocked ejecta ( iLaMassa et al.l 120081 ) . However, 
an additional particle population, older and less energ etic, is needed to explain the mul- 
tiwavelength spectrum of Vela X (Ide Jager et al.l 120081 ) , c omplemented by the recent Ge V 



detections by AGILE (jPellizzoni et al. 



2010|) and Fermi (Ubdo et al.l boioh. The a ngular 



resolution of IBIS/ISGRI on board the INTEGRAL observatory flWinkler et al.l 120031 ) com- 
bined with its large field of view allows for the first time to address the problem of the full 
morphology of the hard X-ray nebula. Here we report on the INTEGRAL identification of 
extended hard X-ray emission from Vela. 
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2. Observations and analysis 
2.1. IBIS/ISGRI imaging 



We analyzed all public IN TEGRAL pointings w ithin 12° from the Vela pulsar. We first 
analyzed the data from IBIS fjUbertini et al.ll2003l ). the coded mask imager on board IN- 
TEG RAL, and in particular of its low energy detector ISGRI (15 keV-1 MeV. iLebrun et al. 
20031 ). The IBIS/ISGRI data have been collected from 1976 pointings between 2003-03 and 
2008-07, for a total exposure time of 5.6 Ms. In the 18-40 keV mo saicked image (Fig. [Q, 
obtained with the Offline Scientific Analysis (IGoldwurm et al.l 120031 ) software v. 8, we found 
a llOcr point-like source at the pulsar position. The point spread function (hereafter PSF) 
encompasses the pulsar, the Chandra PWN, and part of the fainter region. 

An extended emission in the NE/SW direction is also visible in the image, spanning ~50' 
on both sides. After subtraction of the point-like source by fitting it with a 2D Gaussian 
profile (a = 6.2'), the NE side appears larger and more significant than the SW side, which 
is coincident with the ROSAT and H.E.S.S. cocoon (Fig. [TJ second and third panels). The 
extended emission also matches the one found by the Birmingham Spacelab 2 telescope in 
2.5-10 keV (Fig. [[J fourth panel). The individual pixels of the IBIS/ISGRI feature are at the 
~3-6 a significance level. Such a large cluster of low-significance pixels is not observed in the 
rest of the image, and it is not reminiscent of IBIS/ISGRI coding noise. After smoothing, it 
is the only residual excess besides the known point sources. A further evidence of extended 
hard X-ray emission beyond the inner PWN is provided by the spectral analysis. 



2.2. INTEGRAL spectral analysis 



We extracted the IBIS/ISGRI spectrum of the point-like source from mosaicked images 
in narrow energy bands between 18 and 400 keV (Fig. [2]). To account for the evolution of 
the instrument response, we produced an average of the responses weighted by the on-source 
exposures in their r espective vali dity epochs. All the spectra in this work have been fitted 
using Xspec v. 11. 3 f lArnaudlll996l ). and the uncertainties are reported at the 90% confidence 
level. A best fit with a single power law model yielded a photon index Tisgri = 2.00 ± 0.04 
and a flux F = (4.76 ± 0.09) x 10~ n erg cm" 2 s" 1 in the 20-40 keV range (x 2 re d= 1-16/16 
d.o.f.). This spectrum is ~50 times higher than the phase-avera ged one of the Vela pulsar at 
20 keV (Hermsen and Kuiper, priv. comm. ; lHarding et al.ll2002[ ) . Therefore, the IBIS/ISGRI 
emission is dominated by the nebula. 



We also analyzed the data from the INTEGRAL spectrometer SPI (20 keV-8 MeV; 
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Fig. 1. — IBIS/ISGRI significance map in the 18-40 keV range (celestial coordinates, J2000; 
north is up, east is left). The color scale is chosen to highlight the faint emission. First 
panel: significance map; other panels: significance map after subtraction of the point-like 
source, and smoothing with a 3 pixel (<7~7.5') Gaussian kernel. Contours: ROSAT, 0.5- 
2 keV (second panel, green); H.E.S.S., VHE Gamma-rays above 1 TeV (third panel, red , 
Aharonian et al.l 120061 ): Spacelab 2, 2.5-12 keV (fourth panel, cyan, IWillmore et al.lll992l ). 
The cross indicates the pulsar position. The X point in the third panel marks the best-fit 
centre of gravity of the TeV emission. 



Vedrenne et al.l 120031 ) collected simulta neously to the IBIS dat a. SPI spectra have been 
extracted using the SPIROS package Jskinner fc Connelll l2003h within the OSA analysis 
software. The SPI data are best fitted (xl e d = 0.7/10 d.o.f.) by a power law model in 
the 20-300 keV range with photon index Tgpi = 2.15 ± 0.15, compatible within the errors 
with the IBIS/ISGRI one, but with a higher flux, F = (9.1 ± 0.6) x 10~ n erg cm" 2 s" 1 
in 20-40 keV, with respect to IBIS/ISGRI. A joint fit to the IBIS/ISGRI and SPI spectra 
with a single power law yields a compatible photon index, but it requires a renormalization 
constant of 1.8 to recover the higher SPI flux. Such a discrepancy can not be accounted for 
by an intercalibration factor, which is in the 0.8-1.2 range for IBIS/ISGRI and SPI (e.g., 
Jourdain etaJl2008[ iBouchet et al-lbopgh . 



The p hoton index measured by Beppo5y4X/PDS in the same energy range (Ppds = 
2.00 ± 0.05. iMangano et al.ll2005t ) is consistent with both the spectral indices derived above, 
whereas the flux lies between the IBIS/ISGRI and SPI one. As shown in Fig. El the 
IBIS/ISGRI, Beppo&4X/PDS, and SPI fluxes correlate with the respective PSF radii (HWHM: 
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6', 390, and 1.3°), suggesting that each instrument samples a different portion of the neb- 
ula. Due to the coded mask deconvolution of IBIS/ISGRI, optimized for point sources, the 
recon structed flux of an e xtended source of 60' radius is lower than the real one by a factor 
~50 (IRenaud et al.ll2006l ). Therefore, a flux of 4.3 x 10~ n erg cm" 2 s _1 (the difference be- 
tween the SPI and ISGRI fluxes) from such a source would be measured as low as ~10~ 12 
erg cm -2 s _1 by IBIS/ISGRI, close to its sensitivity limit. 

To support the hypothesis that an extended source is present in the IBIS/ISGRI data 
but diluted b y the coded mask de convolution, we refined the analysis following the method 
developed by IRenaud et al.l (120061 ) for analyzing emission from extended sources. 



We extracted the IBIS/ISGRI count rates from concentric circles centered on the pulsar 
with radii up to 80', and converted them into flux by assuming a photon index as for the 
point-like source. This integrated flux as a function of the extraction radius does not reach a 
plateau after 15', as expected for a point-like source, but slowly increases up to ~60' (Fig. [3]). 
The integrated IBIS/ISGRI flux also recovers the BeppoSylX/PDS and SPI fluxes at radii 
comparable with their PSFs. This confirms the detection of extended hard X-ray emission 
beyond the inner PWN. 



2.3. IBIS/ISGRI and Suzaku /XIS combined spectrum 

The MECS and PDS instruments on boa rd BeppoS^X allowed to measure a break 



at energy 12.5 ±1.5 keV (IMangano et al.ll2005l ). However, the different angular resolution 



of the two instruments required the authors to combine the PDS spectrum rescaled by an 
intercalibration factor with the MECS spectrum extracted from a 15' radius region. With 
a much smaller PSF, INTEGRAL IBIS/ISGRI can be combined with Suzaku /XIS on the 
same extraction radius. 

Suzaku observed the Vela pulsar and PWN on 2006 July 10 and 11. Event files from 
version 2.0.6.13 of the Suzaku pipeline were used and spectra were extracted using XSE- 
LECT. Response matrices and ancillary response files were generated for each XIS using 
XISRMFGEN and XISSIMARFGEN version 2007-05-14. The data of the XIS2 camera 
were not considered because of a more uncertain calibration. The effective exposure time 
for each XIS is about 60.3 ks. The source photons are extracted from a circular region with 
a radius of 6' to match the IBIS/ISGRI PSF. Background photons are extracted from blank 



1 As the PDS had no imaging capabilities, this is the HWHM of the instrumental angular response 
(jFrontera et al.lll997t ). 



-6 - 



sky observations within the same region as the source. 



We modeled the pulsar contribution by introducing two black bodies co rrected by the 



inters tellar absorption, with fixed parameters as measured by XMM-Newton (IManzali et al. 



20071 ). An absorbed non-equilibrium plasma emission model (VMEKAL) i s also included 



to acc ount for the thermal supernova remnant; fixing the abundances as in lLaMassa et al. 
( 120081 ). we derived a temperature of k B T = 0.214+q;qq5 keV. The XIS spectra show a bright 



power law compone nt with photon index and flux (Table [TJ compatible with the MECS 
at the same radius ( jMangano et al.l 120051 ) . It is connected to the IBIS/ISGRI spectrum, 
confirming that the IBIS/ISGRI flux is due to the PWN. However, a spectral break is 
required to simultaneously fit the XIS and the IBIS/ISGRI data. The break is located at 
higher energies (27 ± 3 keV) than the one derived by Beppo&4X A simultaneous fit with a 
single power law is statistically rejected (x 2 of 11252 for 6245 d.o.f.). 

We also refined the spectral fit by replacing the broken power law with a Band model, 
an empirical four-parameter model consisting of two power law components smoothly joined 
by an exponential roll-over ( jBand et al.lll993l ). The fitted parameters are compatible with 
the ones derived with the broken power law model (Table [JJ, for a comparable \ 2 ■ Notably, 
the two power laws intersect at 25 ± 7 keV, which corresponds to the break energy in 
the broken power-law model. In the next section we adopt the Band model, as a gradual 
transition should be more representative of the change of slope around a cooling break (e.g., 
KardashevlE962] ). 



3. Discussion 

Thanks to the deep INTEGRAL exposure, we were able to discover diffuse emission 
above 20 keV beyond the inner Vela PWN. Such emission is resolved in two different regions: 
a southern hot-spot coincident with the Vela cocoon, notably with the peak of the TeV 
brightness profile measured by H.E.S.S., and a more extended northern emission, without 
any counterpart at other wavelengths and outside Vela X. However, recent Suzaku/XIS 
observations showed non-ther mal emission below 10 keV in a region located at the boundary 



of the IBIS northern emission (IKatsuda et al.l 1201 if ). We then compared the IBIS spectrum 



to the ones of BeppoSAX /PDS and SPI, and explained the flux differences by their different 
PSFs, which sample different portions of the nebula, confirming a large extension (~1° 
radius) at hard X-rays. 

We also reported on the spectrum within 6' from the pulsar using the IBIS/ISGRI 
and the XIS telescopes. This is the first broad-band X-ray spectrum of the Vela PWN 
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Table 1. Best-fit spectral parameters. 



Broken power law Band 





1 R42 + - 005 

1 - u ^ :Z -0.006 


1 610+ 0012 
i.uiu_ 005 


r 2 


2.07 ± 0.05 


9 9+O.I 
z - z -0.2 


E b (keV) 


27 ± 3 




(keV) 






Ft 


1.11 ± 0.02 


1.11 ± 0.06 


XL (d.o.f) 


1.03 (6235) 


1.03 (6235) 



Note. - Best-fit spectral parameters of the 
combined SuzaAru/XIS-IBIS/ISGRI data within 
6' from the Vela pulsar. The uncertainties are 
at the 90% confidence level. *Eq is the folding 
energy in the Band model. The two power laws 
intersect at (r 2 - 1\) e" 1 E = 25 ± 7 keV, which 
corresponds to the break energy in the broken 
power- law model E b . tFlux in the 20-100 keV 
energy band in units of 10~ 10 erg cm -2 s _1 . 
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taken within a region with the same angular extension below and above 10 keV. The change 
of slope around 25 keV (Ar = 0.59 ± 0.15 for the Band model) is compatible with the 
standard value of 0.5 expected from a cooling break occurring in a continuously injected 
electron distribution affected by radiative losses. Indeed, the XIS and IBIS/ISGRI photon 
indices (~1.6 and ~2.2) are compatible with the synchrotron spectrum of a s hock-accelerate d 



electron distribution in the uncooled and cooled regime, respectively (e.g.. IChevalierll2000l ). 
In this framework, the cooling energy is expected to decrease with time, that is at increasing 
integration radii from the pulsar. The lower break energy measured by BeppoSAX MECS 
and PDS on a larger angular extension may indicate a cooling break propag ating along the 



flow. Indeed, the cooling energy in the cocoon is expected around 1 keV (ILaMassa et al. 



20081 ). The similar photon indices found by IBIS, Beppo&4X/PDS, and SPI suggest that the 
radiative losses above 20 keV already balance the injection rate within the region enclosed 
by the IBIS/ISGRI PSF. 

The measurement of a cooling energy at ~25 keV in the photon spectrum allows us to 
set an upper limit on the time (residence time) spent by the particles in the region within 6' 
from the pulsar, corresponding to a distance of 0.5 <i 2 9o P c - Accounting for synchrotron losses 
and inverse Compton losses in the Thomson regime, the cooling frequency as a function of 
the residence time t can be written 

/ \ SI ml c 9 

Vc{ ) = 327re7 (1 + U ph /U B ) 2 BH*' { ) 

where U p h and Ub are the photon field and magnetic energy densities, respectively. Eq. 
J) coincides w ith the cooling frequency calculated assuming only synchrotron losses (e.g., 



Chevalierl 120001 ) for B ^> ^/9>TiU p h- Solving it for t, the residence time has a maximum 
occurring for a magnetic field B = >/87rU p h/3 independently from the cooling frequency. 
Such a magnetic field amounts to 1.9 fiG when the target photons are provided by the 
Cosmic Microwave Background radiation (CMB, U p h = 0.26 eV cm -3 ). For a cooling energy 
of 25 keV [y c = 6 x 10 18 Hz), the residence time of electrons is 1650 yr in this case, and shorter 
for any different intensity of the magnetic field and any additional photon field. Therefore, 
the electrons radiating in X-rays in the considered region can not be older than 1650 yr. We 
conclude that those injected before have flowed out of the region within 6' from the pulsar. 
This requires a moderate average velocity (> 300<i 2 go km/s). 

We explored the hypothesis that particles of all energies remain in this region for the 
residence time which yields a cooling break at 25 keV for a given B. This is done by means 
of a time-dependent one-zone model of the spectral energy distri bution (S.E.D., see Fig. |H) . 



The injection spectrum is composed by a relativistic Maxwellian (ISironi Sz Spitkovskyll2009f ) 
and a cut-off power law with index 2.2. The radiative losses depend on a magnetic field 
of constant intensity B and a target photon field taken at the Vela position according to 
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Moskalenko et al.l ( 120061 ). The cut-off e nergy is fixed by B through the condition that the 
acceleration rate equals the cooling rate ( jde Jager et al.lll996l ). The shape of the distribution 
is constant, whereas the normalization is proportional to the pulsar spin-down power: 



N 



rjE 



T w m e c 2 ' 



(2) 



where rj is the fraction of E converted into the wind energy, and IV is the average Lorentz 
factor fixed by the energy conservation. Such an injection rate decreases in time following 
E{t) — Eq (1 + t/tdec)' 13 , where E is the initial spin- down power, the s pin-down time- 
scale, and (3 = (n + l)/(n — 1) for a braking index n (IPacini fc Sal vat ill 19 73l ). 



The S.E.D. in Fig. H] is reproduced by this simple model for a range of values of B. 
The upper limit on the integral flux above 1 TeV implies a firm lower limit on the ma gnetic 
field, B > 10 /iG, higher than the one estimated in the cocoon (jde Jager et al.l 120081 ). This 
result depends on the relative intensity of the synchrotron and Compton peaks, which is not 
affected by the conversion efficiency 77 (1.3% for B = 10 /iG). It is also nearly irrespective 
of tdec and n: given the shor t evolution time, a manifold break around the cooling energy 
due to the pulsar spin-down ( Gelfand et al. 2009 ) is not evident^. The extrapolation of the 
X-ray spectrum at lower energies falls squarely on the measured radio fluxes of the compact 
PWN. A fine tuning is obtained adopting a relativistic Maxwellian with a temperature of 3 
GeV for B = 10 /iG; alternatively, a low-energy break i n the injected electron p opulation is 
required at GeV energies, as found in young PWNe (e.g.. lBucciantini et al.ll201ll ). Additional 
measurements between radio and X-ray frequencies are needed to exclude a more complex 
spectrum (e.g.. lSlane et al.ll2008l ). 



The fact that diffuse hard X-ray emission is detected beyond 6' from the pulsar in the 
northern direction strengthens the hypothesis of particle leakage. Unlike the southern hot- 
spot, which may simply be the high-energy counterpart of the X-ray cocoon, the northern 
emission may be gen erated by particles i njected after the transit of the supernova reverse 



shock (~3 kyr ago, iBlondin et al.l 120011 ) . The original PWN is then set apart from the 



pulsar, becoming rel ic, while the latter forms a new PWN in subsonic expansion inside the 
supernova remnant ( Ivan der Swaluw et al.l 120041 ) . To understand the nature of the northern 
emission it will be crucial to extend the X-ray mapping of the Vela PWN with XMM-Newton 
and Suzaku and with the forthcoming focusing telescopes at higher energy X-rays (Nu-STAR 
and ASTRO-H). 



2 We tes ted tdec between 10 yr and 1000 yr, and both n = 3 (dipolar rotator model) and n = 1.6 
(measured. iDodson et al.ll2007l) . Eq has been chosen to yield E — 6.9 x 10 36 erg s^ 1 at the present time. 
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Fig. 2. — Combined Suzaku /XISO-1-3 (blue, green, and red) and IBIS/ISGRI (black) spectra 
of the Vela PWN within 6' from the pulsar fitted to the model with a Band component 
described in Sec. 12.31 The data have been graphically rebinned for clarity. The lower panel 
shows residuals from the best fit in units of la. 
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Fig. 3. — IBIS/ISGRI flux in the 20-100 keV band integrated at varying extraction radii 
(open triangles), along with the BeppoSAX/PDS and SPI fluxes in the same band (filled 
circles). In this case the integration radius is given by their PSF (39', and 1.3° for 
BeppoS^X/PDS, and SPI, respectively). The origin of the x-axis corresponds to the pul- 
sar position. The expected (dotted line) and measured (the bright binary pulsar Vela X-l 
rescaled, open squares) profiles for a truly point source are also shown. 
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Fig. 4. — S.E.D. of the Vela PWN emission within 6' from the pulsar, fitted with the model 
described in the text. The Suzaku XIS and INTEGRAL IBIS/ISGRI spectra are shown in 
red and blue, respective ly. The radio fluxes of the inner Vela PWN are shown (IDodson et al. 
2003I ; lHales et al.l 12004 orange and green circles respectively). The upper limit (99.9%) on 
the integral flux abo ve 1 TeV within 6' from the pulsar is also shown (purple arrow, assuming 
a photon index of 2; lAharonian et al.l 120061 ) . The measurements of th e large-scale PWN are 



reported in gre y for comparison: Vela X in radio (lAlvarez et al 



2001 



Abdo et al.ll2010h . at 



GeV energies ( Abdo et al. 2010 ). and the TeV cocoon ( Aharonian et al. 20061 ) . The total 



(synchrotron and IC) model spectrum is indicated with a thick (thin) solid line. The IC 
emission is computed taking into account the CMB (dashed line), dust (dot-dashed line), 
and star-light (dotted line). The magnetic field is 10 /zG. 



